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ABSTRACT 


A  study  of  the  influence  of  viscosity  ratio  on 
recovery  efficiency  was  carried  out  on  an  eight  feet  long 
unconsolidated  sand  model  employing  fluids  native  to  an 
actual  reservoir  in  the  Lloydminster  "heavy  oil”  area. 

The  average  properties  of  the  two  packs  used  were 
porosity  35«5#j  permeability  3*72  darcies.  Thirteen  tests 
were  run  at  water  saturations  varying  between  13.0#  and 
15*4#.  The  packs  were  preferentially  water  wet. 

Recovery  efficiencies  at  breakthrough  and  WOR  of  5:1 
were  obtained  at  viscosity  ratios  of  292,  646  and  1430. 
Viscosity  ratio  was  controlled  by  temperature  variation  under 
a  back  pressure. 

Rate  studies  were  carried  out  at  viscosity  ratios 
of  292  and  646  which  revealed  a  lowering  of  displacement 
efficiency  at  higher  rates.  This  loss  in  recovery  was 
attributed  to  viscous  fingering.  Breakthrough  recoveries 
at  low  rates  were  found  to  be  29.8#  and  24.6#  of  initial  oil 
in  place  at  viscosity  ratios  of  292  and  646  respectively. 

A  theoretical  recovery  versus  viscosity  ratio  was 
established  using  Buckley -Leverett  theory  and  experimentally 
determined  relative  permeabilities.  Actual  recoveries  were 
higher  than  those  predicted  by  theory. 

Pressure  profiles  were  used  to  observe  the  progress 
of  the  front.  There  was  some  evidence  that  higher  breakthrough 
recoveries  would  be  possible  with  a  longer  model. 
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INTRODUCTION 


The  purpose  of  this  study  was  to  determine  the 

efficiency  of  displacement  of  a  high  viscosity  crude  oil 

from  a  natural  unconsolidated  sand  using  a  formation  water 

as  the  displacing  phase.  The  crude  oil,  water  and  sand 

came  from  wells  in  the  Lloydminster  "heavy  crude"  area. 

This  is  part  of  a  larger  research  project  concerned  with  the 

recovery  of  heavy  oils.  The  impetus  for  a  continuation  of  this 

field  of  research  seems  to  be  assured  as  significant  reserves 

of  high  viscosity  oils  are  found  in  southern  Alberta  and 

southwestern  Saskatchewan  in  addition  to  the  Lloydminster 

area.  It  is  also  possible  that  attempts  may  be  made  in  the 

future  to  produce  the  tar  sands  in  situ.  This  may  be  tried 

in  the  tar  sands  proper  of  the  Athabasca  area'  '  or  in  other 

areas  of  northern  Alberta  where  large  "heavy  crude"  deposits 

(  S'!) 

are  known  to  exist'  -  7 . 

Prior  to  conducting  this  study  a  survey  of  the 
water  flooding  literature  revealed  that  model  experiments 
must  be  properly  scaled  in  order  to  make  the  results  more 
generally  applicable.  Results  of  tests  involving  the 
immiscible  displacement  of  viscous  fluids  suggested  that 
particular  attention  should  be  paid  to  rock  wettability,  rate 
sensitivity,  and  effect  of  viscosity  on  relative  permeabilities. 

Fluids  actually  native  to  a  reservoir  were  used  in 
this  work  in  order  to  simulate  reservoir  sand  wettability. 


Previous  investigations  of  rock  wettability^  °’  66, 

have  shown  it  to  be  a  controlling  variable  in  immiscible  dis- 

(2  3  n  4 ) 

placement.  Other  work'1  *  ’  '  has  shown  that  some  reservoir 

fluids  can  affect  rock  wettability. 

( 13  17  46  SO) 

Scaling  experiments  ’  ’  3  J  '  have  shown  that 

capillary  pressure  can  influence  laboratory  tests  causing 

recoveries  obtained  to  be  rate  sensitive.  However,  by 

increasing  the  injection  rate  a  point  is  reached  above  which 

capillary  effects  are  rendered  negligible.  The  same  effect 

has  been  achieved  at  lower  rates  by  using  longer  systems. 

(  -\  o  i  s  ?6  sp) 

Other  work'  >  >  ■  has  shown  that  recoveries 

may  be  rate  sensitive  at  high  rates  and  high  viscosity  ratios. 
This  rate  influence  has  been  explained  as  being  caused  by 
viscous  fingering. 

A  model,  eight  feet  in  length,  was  designed  in 
order  to  satisfy  the  scaling  requirements  mentioned  above. 

The  role  of  length  in  scaling  considerations  is  discussed 
further  in  a  separate  section  on  theory. 

In  order  to  establish  the  original  wettability  of 
the  sand  some  tests  were  run  on  the  long  core  using  articicial 
liquids  at  a  viscosity  ratio  of  unity.  Relative  perme¬ 
ability  ratios  obtained  from  these  tests  were  used  as  a 
basis  for  determining  the  wettability.  Relative  permeability 
ratios  were  then  obtained  from  tests  using  the  native 
fluids  and  compared . 


-  3  - 

Variation  of  viscosity  ratio  was  achieved  by 
temperature  control.  It  was  known  that  temperature  affects 
interfacial  tension  and  possibly  wettability (3°*  31) .  How¬ 
ever  the  effect  of  changing  interfacial  tension  was  considered 
to  be  minor'  ’  and  any  wettability  change  would  be  reflected 
in  the  unsteady  state  relative  permeability  data^1®). 

In  consideration  of  the  possibility  of  fingering, 
properly  scaled  investigations  in  this  high  viscosity  range 
were  not  available.  For  this  reason  rate  studies  were  carried 
out.  Although  the  exact  nature  of  fingering  was  not  revealed 
in  the  tests,  its  occurrence  was  verified.  Breakthrough 
recoveries  decreased  at  the  higher  rates. 

Steady  state  relative  permeability  curves  were 
obtained  at  the  conclusion  of  the  displacement  tests.  All 
the  relative  permeability  data  were  used  to  construct  a 
theoretical  curve  relating  displacement  efficiency  to 
viscosity  ratio.  The  actual  recoveries  obtained  were  higher 
than  those  indicated  by  the  theoretical  curve.  This  differ¬ 
ence  can  be  explained  by  the  extreme  sensitivity  of  the 
Buckley  Leverett  method  to  small  changes  in  the  relative 
permeability  ratios  in  the  particular  saturation  range  of 
interest . 

It  is  possible  the  scaling  requirements  in  the  model 
were  not  completely  fulfilled  because  of  the  possible  effects 
of  the  temperature  variable.  Until  further  tests  are  run 
this  might  limit  the  direct  application  of  the  results. 
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REVIEW  OF  LITERATURE 

Although  the  main  object  of  this  work  was  to 
determine  the  influence  of  high  viscosity  ratios  on  recovery, 
other  literature  pertaining  to  water  flooding  was  also 
reviewed.  The  work  of  Uren^-^  published  in  1928  is  of 
historical  interest  because  it  discusses  the  effect  of  many  of 
the  variables  that  are  still  being  investigated  today.  These 
include  capillary  forces,  rate,  presence  of  gas,  oil 
viscosity,  interfacial  tension  and  rock  wetting  characteris¬ 
tics.  The  fact  that  the  influence  of  these  parameters  is 
still  not  fully  understood  is  a  testimonial  to  the  complex 
nature  of  fluid  flow  in  porous  media. 

In  1941  Buckley  and  Leverett^?)  developed  a  theory 
which  attempted  to  mathematically  describe  the  immiscible 
displacement  process.  More  recently^)  the  theoretical 
development  has  been  criticized  for  its  lacking  rigor.  How¬ 
ever  its  usefulness  is  at  present  not  questioned. 

The  subject  of  model  scaling  received  notable 
contributions  from  Leverett  in  1941^35)  and  Engleberts  and 
Klinkenberg  in  1949^ -*-7).  Probably  the  most  influential 
paper  on  scaling  problems  was  that  of  Rapoport  and  Leas^-^ 
published  in  1953.  This  paper  concerned  itself  mostly  with 
the  influence  of  capillary  forces  in  laboratory  models.  As 
will  be  discussed  more  fully  in  the  theoretical  presentation 
it  was  found  to  be  advantageous  to  minimize  capillary  effects. 


5 


Most  of  the  laboratory  displacement  tests  since  the  work  of 
Rapoport  and  Leas  have  used  their  scaling  criteria.  Since 
short  models  are  most  convenient  for  laboratory  use,  the 
necessary  scaling  condition  to  attain  minimum  capillary  ef¬ 
fects  has  been  achieved  by  the  control  of  rate. 

A  non -capillary  occurrence  of  equal  importance  in 

model  scaling  is  the  phenomenon  of  viscous  fingering.  Ear  iy 

( 39) 

reference  to  this  condition  was  made  by  Muskat  and 

( 17) 

Engleberts  and  Klinkenberg'  .  A  qualitative  description 
of  fingering  in  a  model  relating  it  to  a  non  uniform  dis¬ 
placement  at  the  inlet  end  of  the  core  was  made  by  Jones- 
Parra  et  al^r'°^.  Further  work  by  Kyte  and  Rapoport 
tended  to  confirm  that  fingering  was  synonymous  with  this 
inlet  end  effect" .  This  effect  was  explained  as  being 
localized  imbibition  of  water  at  the  inlet  core  face.  It 
was  suggested  that  fingering  was  aggravated  by  short  systems, 
high  water  injection  rates  and  high  oil  to  water  viscosity  ratios. 
A  quantitative  analysis  of  displacement  from  a 

different  approach  was  made  by  Chuoke  et  al  in  195^ . ^  . 

Fingering  was  created  at  displacement  interfaces  between 
parallel  plates  in  a  form  which  could  be  theoretically  pre¬ 
dicated.  This  theory  was  also  applied  with  reasonable  success 

to  oil  displacement  in  a  medium  of  glass  beads.  Later 

(  ]  7 ) 

work  by  de  Haanv  '  substantiated  the  theory  and  suggested 
that  it  applied  as  well  to  water  wet  systems  with  an 
initial  connate  water  saturation.  Experimental  evidence 


*5*. 
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suggested  that  fingering  would  not  occur  in  oil  wet  systems. 
Later  theoretical  developments  by  0utmans^5)  outlined 
revisions  to  Chuoke's  theory  and  suggested  that  viscous 
fingering  in  the  reservoir  would  have  a  negligible  effect. 

There  are  a  few  papers  which  were  concerned  with 
displacement  efficiencies  at  high  viscosity  ratios  i.e. 
greater  than  one  hundred  to  one.  Of  these,  the  existence  of 
fingering  in  some  form  was  recognized  by  Kyte  and  Rapoport^2)  ^ 
Newcombe  et  al^1^  and  Croes  and  Schwarz^2).  An  interesting 
illustration  of  the  effect  of  fully  developed  fingering 
can  be  observed  by  comparing  the  latter  work  with  that  of 
de  Haan.  This  was  not  pointed  out  by  de  Haan  but  it  seems 
to  be  a  valid  conclusion  as  the  same  model  was  used  in  both 
studies.  Part  of  the  work  of  Newcombe  et  al  showed  the 
effects  of  rate,  interfacial  tension  and  contact  angle  at  a 
high  viscosity  ratio.  They  found  that  high  rates  caused 
decreased  recovery  at  breakthrough  in  a  water  wet  system 
and  that  a  reduction  of  interfacial  tension  gave  slightly  lower 
breakthrough  recoveries  but  higher  ultimate  recoveries . 

Two  previous  works  by  Felsenthal^ 19)  and  Fried^20) 
were  concerned  with  the  influence  of  viscosity  ratio  on 
recovery.  Felsenthal  did  not  deal  with  viscosity  ratios 
above  one  hundred  to  one.  Fried  worked  with  a  short  alundum 
core  and  did  not  consider  the  possible  occurrence  of  fingering. 

Moore  and  Slobod^S)  aid  extensive  work  on  the 
effect  of  rate,  wettability  and  core  length  on  recovery. 


« 
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They  concluded  that  water  wet  cores  did  not  offer  any 
particular  scaling  problems.  They  observed  that  recovery 
histories  of  one  inch  core  plugs  were  similar  to  those  of 
long  cores  when  flooded  under  the  same  conditions  of  rate, 
viscosity  level  and  interfacial  tension.  It  may  be  signifi¬ 
cant  that  though  they  varied  the  viscosity  of  the  displaced 
fluid  up  to  144  centipoises,  the  viscosity  ratio  was  kept 
at  approximately  one. 

It  has  been  a  fairly  common  practice  in  relative 
permeability  studies  to  use  viscous  oils  for  displacement 
tests  in  order  to  extend  unsteady  state  relative  permea¬ 
bility  data  to  a  larger  saturation  range.  Referring  to 
this  practice,  Loomis  and  Crowell^^)  concluded  that  the 
displacement  method  was  not  entirely  trustworthy  with  water 
wet  cores  in  the  absence  of  check  experiments. 

Odeh(42)  formulated  a  theory  which  suggested  that 
relative  permeability  was  a  function  of  viscosity  ratio. 
Experimental  data  obtained  with  fairly  high  viscosities 
included  non  wetting  phase  relative  permeabilities  of  values 
much  greater  than  one.  Odeh  suggested  that  media  of  permea¬ 
bility  greater  than  one  darcy  would  not  exhibit  this  phenomenon. 
Considerable  controversy  has  arisen  since  this  theory  was 
formulated ( ^  . 

Much  of  the  previous  work  on  models  has  involved 
the  use  of  artificial  fluids  to  avoid  the  possible  influence 
of  impurities^^  and  polar  compounds^2*  .  There  has  been  a 


p  „ 


.recent  trend  toward  simulating  reservoir  conditions  as 
closely  as  possible^’  in  order  to  reproduce  the 

proper  conditions  of  wettability.  It  was  observed  that  in 
many  cases  cores  which  exhibited  oil  wet  tendencies  at 
room  conditions  would  become  more  water  wet  under  the  in¬ 
fluence  of  increased  temperature  and  pressure. 


-  9  - 

THEORY 


SCALING  OF  CAPILLARY  FORCES 

The  fractional  flow  formula  was  first  developed 
by  Buckley  and  Leverett^?)  for  an  immiscible  displacement.  If 
water  is  displacing  oil  in  a  horizontal  system  the  equation 


becomes 


1  - 


- 


K0  (dPc) 

Mo  Qt 


1  j  ^o  /Uw 
kw  Ao 


(1) 


where:  f  ,  „ 

w 


ko*  - 


qt - 


L  - 


fraction  of  water  flowing  at  any  point 
effective  permeability  to  oil  (darcies) 
relative  permeability  to  oil  and  water 
respectively 

viscosity  of  oil  and  water  respectively 
(centipoise) 

total  flow  rate  (cc/sec) 
capillary  pressure  (atm) 
length  (cm) 


The  influence  of  capillary  forces  can  be  realized 

by  studying  the  right  hand  term  in  the  numerator  of  equation 

(l).  It  can  be  shown  that  regardless  of  the  wettability  of 

dPc 

the  system  the  capillary  pressure  gradient  —  is  always 

dL 

negative.  This  produces  a  positive  contribution  to  the 
numerator  indicating  a  larger  fraction  of  water  flowing  at 


+ 


1.0 


points  where  there  is  a  capillary  pressure  gradient.  One 
would  expect  that  an  increase  in  oil  viscosity  or  rate  should 
decrease  the  capillary  effect  because  they  appear  in  the 
denominator  of  equation  (l). 


/  r-j  \ 

These  conclusions  were  confirmed  by  Rapoport'  ' 
who  developed  flow  equations  in  dimensionless  form  to 
arrive  at  a  "scaling  coefficient".  It  was  found  that  a 
critical  scaling  coefficient  could  be  established  for  a 
given  system  above  which  capillary  forces  would  not  decrease 
the  breakthrough  recovery  efficiency.  In  a  later  paper^ 50) 
it  was  verified  that  increasing  the  oil  viscosity  would 
decrease  the  value  of  the  critical  scaling  coefficient  other 
quantities  remaining  the  same. 

(  32) 

From  the  evidence  of  Kyte  and  Rapoport'  '  the 
wetting  phase  saturation  build  up  at  the  outlet  core  face 
(sometimes  called  outlet  end  effect)  is  negligible  provided 
that  capillary  forces  are  scaled  in  the  model. 

The  dimensionless  capillary  scaling  number  I  as 
used  by  de  Haan  who  assumed  a  zero  contact  angle  ist 


where:  L 


l  y~Mu 
a-i/T1 


length  (cm) 


w 


water  viscosity  (poise) 


~V~  -  superficial  velocity  (cm/sec) 
T  -  interfacial  tension  (dynes/cm) 

p 

K  -  permeability  (cm  ) 
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The  numerator  in  this  expression  is  Rapoport’s 
scaling  coefficient  when  proper  regard  is  given  to  the 
choice  of  units . 


VISCOUS  FINGER  SCALING 


Chuoke’s  scaling  theory  is  useful  to  this  work 
despite  the  theoretical  objections  of  Outmans^^)  because 
it  is  adaptible  to  a  semi  empirical  approach. 

The  following  relation  was  derived  by  Chuoke 
et  al^1(^  for  the  average  distance  between  viscous  fingers. 

a  _  c 

/m  -  hu/0  -A*) 

where:  average  wave  length  (cm) 

^T"  interfacial  tension  (dynes/cm) 

K  permeability  (cm2) 

V*  superficial  velocity  (cm/sec) 
and/(y  oil  and  water  viscosity  (poise) 

C  constant  (varies  with  system) 


The  smallest  wave  length  possible  at  the  onset  of 
fingering  is  called  the  critical  wave  length  ^  cr  which  is 
related  to  the  average  distance  between  fingers  by 


^cr  =  Am//T"  W 

No  fingering  will  exist  when  $cr  is  greater  than 
the  largest  lateral  dimension  of  the  system.  In  applying 
this  criterion  to  a  reservoir,  de  Haan  used  the  thickness 
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of  the  reservoir  as  the  reference  dimension  of  the  system. 
This  seems  to  be  a  point  of  possible  controversy  but  has 
been  accepted  for  this  development.  The  displacement 
efficiency  at  breakthrough  is  a  function  of  the  number  of 
fingers  which  is  in  turn  a  function  of  the  rate  of  displace*- 
ment.  When  about  five  or  more  fingers  have  formed  the 
system  is  no  longer  rate  sensitive,  and  recovery  efficiency 
at  breakthrough  stabilizes  at  some  lower  level.  Fingering 
does  not  appear  to  occur  in  oil  wet  systems. 

The  constant  C  takes  into  account  the  fact  that 
the  interfacial  tension  at  the  interface  in  a  porous  medium 
containing  a  connate  water  will  be  something  greater  than 
that  measured  at  a  planar  interface,  de  Haan  suggests  a 
C  value  of  about  300  when  connate  water  is  present. 

If  the  onset  of  fingering  is  observed  in  a  given 
system,  C  can  be  evaluated  because  at  this  point  ^cr  must  be 
equal  to  h  so  that 

C  r 

where  in  the  case  of  a  core,  h  is  the  diameter  in  cms. 

COMPATIBILITY  OF  CAPILLARY  NUMBER  AND  FINGER  SCALING 

The  application  of  the  capillary  number  and  finger 
scaling  to  a  model  might  appear  straight  forward.  It  is 
required  to  make  these  quantities  the  same  in  the  model  as 
in  the  reservoir  without  violating  other  scaling  groups. 


ith  ci 

k 


(5) 
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The  most  convenient  parameter  to  adjust  in  the  model  is 

rate  and  this  can  be  best  illustrated  with  the  use  of  examples. 

Case  Is  Reservoir  conditions:  Fingering  absent 

Stabilized  zone  negligible. 

It  is  necessary  to  find  a  rate  for  the  model  which 
would  result  in  a  capillary  number  greater  than  some  minimum 
value,  de  Haan  suggests  1^0. 1.  If  this  rate  causes 
fingering  it  would  be  impossible  to  properly  scale  the  test 
without  re-designing  the  system. 

The  point  to  be  made  here,  which  is  applicable  to 
this  work,  is  the  desirability  of  having  a  long  core  when 
finger  scaling  is  an  unknown  factor;  i.e.  C  is  unknown. 

The  model  length  does  not  enter  finger  scaling  directly  but 
a  long  core  permits  capillary  scaling  criteria  to  be  fulfilled 
at  lower  rates. 

Case  II:  Reservoir  conditions:  Fingering  present 

Stabilized  zone  negligible. 

This  does  not  involve  a  scaling  problem  because  a 
rate  can  likely  be  set  to  scale  the  fingers  and  it  is  reasonable 
to  assume  that  at  this  rate  capillary  effects  will  be  absent. 

However  in  order  to  hypothesize  the  existence  of 
fingering  in  the  reservoir  a  C  value  must  be  available  from 
some  other  experimental  data.  Therefore  it  might  be  advisable 
to  carry  out  a  rate  study  to  verify  this  value. 
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EX  PER I KENT A L  EQU I PMENT 

The  core  holder  was  an  8-foot  section  of  schedule 
80  steel  pipe  of  1-1/2  inches  inside  diameter.  A  1-1/2  inch 
series  900  flange  was  welded  to  each  end.  The  flanges  were 
made  up  using  a  brass  gasket  dressed  with  Garlok  Sealing 
Compound  No.  101. 

The  sand  was  held  in  place  by  sintered  bronze 
discs  which  were  tapped  snugly  into  place,  flush  with  the 
end  of  the  core.  The  presence  of  the  flange  gasket  insured 
a  void  space  sufficient  for  good  communication  across  the 
core  face. 

Pressure  points  were  located  every  foot  along  one 
side  of  the  pipe.  Small  sintered  bronze  discs  located  above 
the  1  8  inch  ports  prevented  loss  of  sand  from  the  pipe. 

Two  different  views  of  the  entire  apparatus  can  be 
seen  in  Figure  1. 

The  entire  core  assembly  was  set  inside  a  rectangular 
metal  box  long  enough  to  accomodate  inlet  and  outlet  con¬ 
nections  and  deep  enough  to  completely  submerge  the  flanges 
in  a  bath  oil.  The  oil  was  circulated  continuously  and  the 
temperature  was  controlled  to  within  ztl.0°F.  The  primary 
heat  source  was  an  electrical  element.  Additional  piping 
within  the  bath  was  used  to  circulate  water  as  necessary 
to  expedite  heating  or  cooling. 


* 

»  / 
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Figure  1:  Experimental  Equipment 
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A  model  2249  W  II  double  barrelled  constant  rate 
Ruska  pump  was  used  for  all  tests.  Shown  In  the  top  photo  of 
Figure  1,  are  the  stainless  steel  reservoirs  which  were  used 
to  hold  the  brine  and  dead  crude. 

The  collecting  graduates  were  kept  in  a  settling 
bath  controlled  at  the  same  temperature  as  that  of  the  core. 

A  schematic  diagram  of  the  equipment  is  shown  in 
Figure  2.  Provision  for  the  use  of  live  crude  was  made 
although  this  was  not  used  for  the  tests.  Varsol  was  used  in 
the  lines  to  transmit  pressure  to  the  Heise  gauges. 


AIR  SUPPLY 
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EXPERIMENTAL  PROCEDURE 

SOURCE,  STORAGE  AND  HANDLING  OP  MATERIALS 

Natural  sand  from  a  lease  storage  tank  was  used. 

The  sand  was  first  washed  with  hot  tap  water  and  then 
flushed  with  Varsol.  The  sand  was  then  given  alternate 
washes  of  benzene  and  distilled  water  until  all  traces  of 
staining  were  removed.  Some  fines  may  have  been  lost  during 
these  washings.  The  sand  was  then  placed  in  a  drying  oven 
for  a  few  hours  at  about  250°F. 

The  crude  oil  used  in  all  tests  was  obtained  after 
treatment  for  removal  of  water  only.  Crude  oil  viscosity  was 
measured  at  room  conditions  with  a  Brookfield  rotating  bob 
type  viscometer  using  as  calibrating  oils  SAE  20,  30  and  50* 
Viscosities  at  higher  temperatures  and  pressures 
were  obtained  from  tests  on  a  14  foot  piece  of  1/16  inch 
inside  diameter  copper  tubing  immersed  in  the  temperature 
bath.  A  constant  for  this  system  was  obtained  based  on 
the  crude  viscosity  at  room  conditions.  This  constant  was 
used  to  convert  pressure  drops  at  other  conditions  into  vis¬ 
cosities.  Welker  and  Dunlop(57)  describe  a  more  elaborate 
set  up  for  this  type  of  equipment. 

The  formation  water  was  filtered  through  filter 
paper  and  stored  in  a  glass  container.  Considerable  mist 
was  removed.  It  was  usually  necessary  to  filter  the  brine 
again  when  it  was  transferred  to  the  lucite  storage  reservoir. 
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The  brine  was  de-aerated  before  using. 

Surface  and  interfacial  tension  measurements 
were  made  with  a  Du  Nouy  ring  tensiometer.  As  the  constant 
temperature  bath  was  separate  from  the  instrument  the  con¬ 
tainer  was  brought  to  bath  temperature  and  quickly  moved  to 
the  instrument  platform.  The  temperature  of  measurement 
was  taken  as  an  average  of  the  bath  temperature  and  the 
temperature  obtained  immediately  after  the  ring  pulled  free. 

The  measured  fluid  properties  are  shown  in  Figures 
1-1  to  1-5  and  table  1-1  of  Appendix  I.  Attention  is  drawn 
to  the  fact  that  the  injected  and  effluent  fluids  had  dif¬ 
ferent  interfacial  tensions  presumably  due  to  the  filtering 
action  of  the  sand  pack. 

PACKING  THE  CORE 

One  end  of  the  core  holder  was  completely  flanged 
up  with  the  bronze  screen  in  place  and  a  bull  plug  inserted 
in  the  port.  The  core  holder  was  placed  in  a  vertical 
position  and  half  filled  with  distilled  water.  While  dry 
sand  was  poured  into  the  top  using  a  funnel,  the  base  of  the 
pipe  was  vibrated  using  a  1/2  H.P.  electric  hammer,  Within 
an  hour  the  water  level  rose  to  the  level  of  the  upper  flange. 

At  this  point  a  piece  of  1-1  2  inch  outside  diameter  lucite 
about  1-1/2  feet  long  was  inserted  inside  the  top  of  the  core 
holder.  The  sand  was  then  poured  through  the  lucite  and  when 
the  sand  level  did  not  change  with  continued  vibration  throughout 
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the  entire  length  of  the  core  holder  the  packing  was 
terminated . 

DETERMINATION  OF  FUNDAMENTAL  CORE  PROPERTIES 

To  obtain  the  porosity  the  core  was  dried  for  one 
day  at  a  temperature  of  130°F .  Operating  the  Ruska  pump 
by  hand,  de-aerated  brine  was  injected  into  the  evacuated 
core  until  fully  saturated.  The  pore  volume  was  then  cal¬ 
culated  from  the  pump  readings. 

Prior  to  brine  saturation  of  the  core  an  air 

permeability  was  run  at  room  temperature.  The  upstream  and 

downstream  pressure  taps  were  not  located  at  the  end 

( ?1^ 

faces  as  recommended  by  Hamilton v  ’  .  The  effect  of  this 
on  the  obtained  pressure  profile  could  not  be  ascertained 
because  of  the  variation  in  permeability  inherent  in  the 
packing  procedure.  Following  saturation  of  the  core  with 
brine,  permeability  profiles  at  100$  liquid  saturation  were 
recorded.  In  the  case  of  pack  number  5  many  of  these 
profiles  were  run  at  various  times  using  Varsol  as  well  as 
brine  as  the  saturating  liquid.  The  sand  and  test  model 
properties  are  shown  in  Figure  3. 

Initial  Imbibition  tests  were  run  on  small  sand 
cells,  saturated  with  Varsol  and  immersed  in  distilled  water. 
The  results  indicated  the  sand  was  water  wet.  When  Lloyd- 
minster  crude  was  used  as  the  saturating  oil  no  imbibition 
was  observed.  This  was  attributed  to  the  effect  of  viscosity 


on  imbibition  rate. 
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Five  displacement  tests  were  run  on  core  pack 
number  5  using  Varsol  and  distilled  water.  The  objectives 
sought  were  to  study  the  relative  displacement  efficiencies 
of  oil  and  water  as  the  displacing  phase,  and  to  obtain  an 
unsteady  state  relative  permeability  curve  which  could  be 
used  as  a  reference  for  similar  curves  obtained  from  later 
displacement  data.  In  two  of  the  five  tests  the  water  was 
treated  with  sodium  carbonate  to  promote  a  water  wet 

condition^2^) , 

Steady  state  relative  permeability  tests  were  run 
at  the  conclusion  of  the  displacement  tests.  The  range  of 
water  oil  ratios  selected  varied  from  1:3  to  15:1.  The 
total  flow  rate  was  always  close  to  30  cc's  per  hour.  At 
the  conclusion  of  the  last  ratio  the  rate  was  dropped  to 
10  cc*s  per  hour  to  check  for  possible  rate  sensitivity.  No 
provision  was  made  for  differentiating  between  the  pressure 
drop  in  the  water  and  crude  phases.  Saturation  equilibrium 
within  the  core  was  considered  attained  when  injected  and 
effluent  water  oil  ratios  were  identical  and  the  pressure 
profile  was  reasonably  constant.  The  steady  state  relative 
permeability  data  are  shown  in  Tables  II-l  to  II -4  of 
Appendix  II  along  with  a  sample  calculation. 

DISPLACEMENT  TESTS 

Sixteen  displacement  tests  were  carried  out  with 
the  core  saturated  at  some  initial  brine  saturation.  These 
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tests  are  numbered  6  to  21.  Tests  6  to  8  involving  an  11 $ 
water  emulsified  crude  do  not  form  a  part  of  this  thesis. 

Tests  9  to  13  inclusive  were  run  to  study  the  effect  of 
rate  at  the  viscosity  ratio  obtained  at  100°F.  Conditions 
between  tests  9  and  10  were  not  changed.  The  core  was 
repacked  after  Test  13.  Tests  14  to  1 6  were  carried  out  on 
pack  number  6  at  a  constant  rate  of  30  cc's  per  hour  and 
the  viscosity  ratios  obtained  at  temperatures  of  77°F,  100°F 
and  130 °F .  Test  17  to  21  comprised  a  rate  study  at  130°F. 

The  direction  of  flow  was  reversed  for  Test  21. 

When  saturating  the  core  with  crude  it  was  found 
that  very  little  water  was  displaced  after  breakthrough  of 
the  crude  at  the  outlet  end  and  an  irreducible  water  satura¬ 
tion  no  greater  than  13$  was  indicated  for  the  sand.  During 
a  displacement  the  average  water  saturation  of  the  core 
would  rise  to  about  40$  or  greater  depending  on  the 
cumulative  brine  injected.  This  value  could  be  restored 
to  its  initial  value  by  the  injection  of  about  one  pore 
volume  of  crude  at  the  conclusion  of  a  displacement.  This 
procedure  for  resaturation  has  been  recommended  by  Maguss(37) . 

While  reinjecting  crude  after  Test  number  15  the 
bath  temperature  was  increased  from  77°P  to  130°F.  The  volume 
of  effluent  fluid  exceeded  the  injected  volume  by  25  cc.  This 
amount  was  in  line  with  that  expected  from  thermal  expansion. 

A  pressure  profile  was  recorded  prior  to  each  test 
to  confirm  re-establishment  of  saturation  equilibrium  within 
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the  core.  While  obtaining  a  profile  during  the  displacement 
the  front  would  move  a  significant  distance  depending  on 
the  rate.  Therefore  a  time  interval  was  recorded  for 
each  profile  and,  in  addition,  the  exact  time  the  pressure 
was  recorded  at  the  tap  most  recently  passed  by  the  front. 
Profiles  obtained  at  high  rates  had  fewer  anomalous  points 
than  those  at  low  rates.  At  the  lowest  rate  of  2.5  cc’s  per 
hour  the  front  was  difficult  to  discern  from  the  pressure 
profile.  This  was  attributed  to  the  fact  that  the  capillary 
pressure  was  of  the  order  of  the  applied  pressure  gradient. 
Pressure  data  are  shown  in  Tables  III-l  to  IIX-13  of 
Appendix  III. 

Breakthrough  was  anticipated  fairly  well  by  observing 
the  pressure  profiles.  Just  prior  to  the  appearance  of 
large  discrete  water  drops  in  the  effluent,  the  downstream 
pressure  would  fall  slightly  due  to  an  effect  at  the  back 
pressure  regulator.  The  most  diagnostic  evidence  was 
obtained  using  a  sonic  viscometer  with  a  vibrating  probe. 

For  some  reason,  possibly  emulsification  brought  about  by 
the  high  shear  rate,  the  first  appearance  of  water  was 
signalled  by  a  sharp  increase  in  the  viscosity  reading. 

This  instrument  was  used  for  Tests  9  to  13.  The  use  of 
this  instrument  entailed  the  introduction  of  a  cell  of  7 
cc  volume  at  the  end  of  the  core.  Because  of  the  channeling 
action  of  the  water  past  the  crude  in  the  cell,  corrections 
to  the  recoveries  were  difficult  to  administer  for  which  no 
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satisfactory  solution  was  found.  For  this  reason  this 
instrument  was  not  used  for  Tests  14  to  21. 

The  water  in  the  effluent  was  separated  from  the 
oil  in  a  settling  bath.  Recoveries  with  small  water  cuts 
were  centrifuged  to  make  sure  that  volumes  measured  after 
settling  alone  were  not  in  error  due  to  emulsification.  There 
was  a  tendency  for  droplets  of  crude  to  cling  to  the  sides 
of  the  graduates  if  they  were  not  thoroughly  cleaned.  There¬ 
fore  glassware  was  rinsed  with  concentrated  hydrochloric 
acid  when  cleaning.  Recoveries  for  the  tests  are  shown  in 
Tables  III-14  to  III-26  of  Appendix  III. 

To  clean  the  core  the  temperature  was  raised  to 
130°F  to  reduce  the  crude  viscosity.  Using  air  pressure, 
Varsol  was  first  injected  through  the  pressure  taps  until 
free  flow  was  obtained.  Varsol  then  was  injected  at  the 
upstream  end  until  the  effluent  was  clear.  This  required 
about  5  gallons  of  Varsol.  This  was  followed  by  about  a 
gallon  of  pentane,  which  was  then  blown  out  with  dry  air. 

INTERNAL  CORROSION 

When  repacking  the  core  the  bronze  end  screens  as 
well  as  the  inside  of  the  core  holder  were  observed  to  be 
covered  by  a  black  deposit.  This  appeared  to  be  a  corrosion 
product . 
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DISCUSSION  OP  RESULTS 


MODEL  PROPERTIES 

Referring  to  Figure  3  it  will  be  noticed  that 
although  pack  number  6  had  a  slightly  coarser  grain  size 
distribution  than  pack  number  5  the  latter  had  a  slightly 
higher  porosity  and  permeability.  This  may  be  explained 
by  a  difference  in  the  packing  technique  employed  in  each 
case.  Pack  number  6  was  packed  in  a  period  of  about  3  hours. 
The  packing  of  sample  number  5  was  interrupted  for  a  period  of 
days  and  a  smaller  electric  hammer  was  used.  This  also  may 
explain  why  pack  number  5  had  a  minimum  permeability  near 
the  middle  instead  of  at  the  top  as  was  the  case  of  pack 
number  6. 

PRELIMINARY  TESTS 

The  recovery  histories  for  the  five  displacements 
at  unity  viscosity  ratio  are  shown  in  Figure  4.  Higher 
breakthrough  recoveries  were  obtained  when  displacing  Varsol 
with  water.  The  recovery  histories  of  the  treated  and 
untreated  water  displacements  were  reasonably  similar.  The 
sodium  carbonate  treatment  was  intended  to  insure  a  water 
wet  sand^2^  . 

Relative  permeability  ratios  obtained  from  these 
tests  are  shown  in  Figure  5*  These  were  calculated  using 
the  Welge  modification  to  the  Buckley -Leverett  theory.  Some 
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points  obtained  subsequently  from  tests  using  crude  and 
native  brine  are  shown  as  well.  It  was  possible  to  establish 
a  representative  line  through  these  data  points.  Also  shown 
on  this  figure  is  an  area  covering  the  wettability  spectrum 
as  established  by  Patt  and  Klikoff^1®)  for  a  pack  of 
reasonably  similar  grain  size  distribution.  A  comparison 
of  the  curves  would  suggest  that  pack  numbers  5  and  6  were 
water  wet.  The  data  points  of  the  preliminary  tests  were 
obtained  from  a  drainage  direction  saturation  change.  It 
was  necessary  to  use  these  points  to  establish  a  curve  in 
the  saturation  range  of  interest. 

The  recovery  and  pressure  histories  for  two  of  the 
preliminary  tests  are  compared  with  recovery  and  pressure 
histories  for  two  different  wetting  conditions  as  obtained 
from  the  work  of  de  Haan^1^)  in  Figure  6.  The  "water 
displacing  Varsol"  curve  of  this  work  resembled  the  inter¬ 
mediate  wettability  case  in  de  Haan’s  work.  Intermediate 
wettability  was  defined  by  de  Haan  as  being  weakly  prefer¬ 
entially  water  wet. 

DISPLACEMENT  TESTS 

The  discussion  of  the  tests  involving  the  crude 
oil  of  high  viscosity  will  be  concerned  with  uniformity 
of  initial  conditions,  actual  recoveries  obtained,  evidence 
of  fingering,  ideas  on  the  subordinate  production  mechanism, 
theoretical  recovery  and  comparison  of  the  theoretical 
recoveries  with  those  obtained  from  the  literature. 


BRINE  PERMEABILITY  PROFILES 
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Initial  Conditions.  A  summary  of  the  initial  test  conditions 
is  given  in  Table  I.  It  will  be  noted  that  initial  water 
saturation  was  controlled  within  reasonable  limits  and  the 
small  variations  did  not  appear  to  correlate  with  any  other 
factor.  This  was  not  true  of  the  values  for  the  relative 
permeability  to  oil  at  connate  water  saturation  however. 

There  are  three  fairly  obvious  trends  perceivable.  Firstly 
pack  number  5  exhibited  higher  effective  oil  permeabilities 
than  pack  number  6.  Secondly  the  effective  permeability 
decreased  as  the  temperature  increased.  Since  most  of  the 
runs  at  low  temperature  were  made  on  pack  number  5  the 
differences  might  have  been  due  to  either  packing,  or  tempera¬ 
ture  or  both.  Thirdly,  the  effective  permeability  to  oil 
increased  with  increased  rate. 

These  variations  can  possibly  be  explained  by  the 
Jamin  effect(3^)t  At  higher  pressure  drops  the  connate 
water  held  by  capillary  forces  was  forced  to  move,  increasing 
the  path  available  for  flow.  This  is  substantiated  to  some 
extent  by  the  fact  that,  water  cuts  of  about  2 %  were  obtained 
in  the  primary  production  at  higher  rates  and  only  a  trace  of 
water  at  lower  rates. 

The  theory  of  Odeh^2)  would  not  explain  the 
relative  permeabilities  greater  than  one.  Odeh  suggested 
1  Darcy  as  an  upper  limit  to  the  permeability  range  in  which 
viscosity  influences  relative  permeability.  Odeh's  actual 
tests  were  made  with  viscosities  of  less  than  100. 
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In  an  effort  to  explain  the  anomalous  effective 
oil  permeability,  profiles  were  plotted  for  the  tests  at 
different  temperatures.  These  are  shown  in  Figure  7.  While 
there  was  an  apparent  temperature  effect  between  100°F  and 
130°F  the  profile  at  77°F  did  not  confirm  this. 

Figure  8  shows  two  profiles  obtained  for  identical 
conditions  except  direction  of  flow.  These  profiles  were 
run  primarily  to  see  if  there  was  any  kind  of  saturation 
gradient  inherent  in  the  resaturation  technique  employed. 

The  fact  that  the  core  was  removed  from  the  bath  to  turn  it 
end  for  end  and  thereby  allowed  to  cool  and  reheat  may  have 
in  some  way  caused  the  very  large  shifts  evident. 

Recovery .  Recoveries  are  summarized  in  Table  I.  The  recovery 
histories  are  shown  in  Figures  III-3,  III-4  and  III-5  of 
Appendix  III.  Test  numbers  9  and  10  were  run  under  the  same 
conditions  and  give  assurance  of  reasonably  good  reproduci¬ 
bility.  Test  number  14  run  on  pack  number  6  gave  recoveries 
in  line  with  those  of  the  previous  pack  when  consideration 
was  given  to  rate. 

Test  number  21  showed  a  lower  breakthrough  recovery 
than  test  number  20  run  under  identical  conditions  except  for 
the  direction  of  displacement.  This  may  have  been  because  the 
rate  used  for  these  tests  was  in  the  range  where  fingering 
was  apparent.  The  good  agreement  between  tests  for  the 
recoveries  at  a  WOR  of  5:1  suggested  that  the  poor  reproduci¬ 
bility  of  breakthrough  recoveries  might  have  been  due  to  a  rate 
Induced  instability. 
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Figure  9  shows  the  recovery  efficiency  as  a  function 
of  the  capillary  scaling  number.  The  injection  rate  is  also 
shown  in  an  approximate  position  along  the  horizontal  axis.  It 
may  be  observed  that  there  is  a  decline  in  recovery  with  rate 
for  the  viscosity  ratios  of  292  and  646.  For  the  viscosity 
ratio  of  1430  only  one  rate  was  used. 

Upon  an  examination  of  these  results,  two  observations 
were  made.  Firstly,  at  low  rates  the  recovery  approached  a 
constant  value.  Secondly  the  scaling  number  I  at  the  lowest 
rate  employed  was  approximately  0.05.  This  is  less  than  the 
critical  scaling  number  of  0.1  suggested  by  de  Haan. 

A  point  has  been  selected  on  each  curve  of  the 
figure  to  represent  the  rate  at  which  fingering  has  started 
to  affect  recovery.  The  selection  of  this  point  is  very 
arbitrary.  Chuoke’s  C  value  which  corresponds  to  these 
selected  points  is  250  as  indicated  in  Table  .2  To  give  an 
idea  as  to  the  range  of  values  which  C  could  take,  values 
have  been  calculated  for  all  rates. 

Figure  10  compares  the  recovery  history  of  test  number 
9  of  this  work  with  that  of  a  similar  model  from  reference  (5)* 
There  is  a  small  difference  in  viscosity  ratios  between  these 
tests.  A  higher  breakthrough  recovery  was  obtained  in  this 
work  but  further  into  subordinate  production  there  is  reasonable 
agreement.  Test  number  9  was  affected  by  fingering.  A  compari¬ 
son  can  also  be  made  of  the  recoveries  obtained  in  this  work  at 
a  viscosity  ratio  of  292  with  that  of  Newcombe  et  al^1^  . 
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FIG-  10  COMPARISON  OF  RECOVERY  DATA  FOR  TEST  NO- 9  OF  THIS 
WORK  WITH  UNCONS-  SAND  MODEL  QF  CARPENTER  ET  AL  15) 


Evidence  of  Fingering.  The  recovery  curves  shown  in  Figure  9 
are  similar  to  those  of  de  Haan^1^)  except  that  in  the  work  of 
the  latter  a  stable  region  of  fully  developed  fingering  was 
reached.  If  we  assume  that  stability  is  attained  at  an  h/PVm 
value  of  five  or  six,  a  rate  twenty-five  to  thirty-six  times 
the  critical  would  be  required  to  achieve  fully  developed 
fingering.  This  rate  would  have  exceeded  the  pressure  for 
which  the  core  holder  was  designed.  Had  this  been  possible 
it  might  have  given  some  support  for  Chuoke's  theory  as 
opposed  to  some  other  effect. 

In  an  effort  to  further  analyze  the  fingering 
mechanism  the  data  were  plotted  to  provide  fractional  flow 
curves,  pressure  profiles,  injectivity  histories  etc. 

The  fractional  flow  curves  obtained  are  shown  in 
Figures  11  and  III-l.  If  there  were  no  rate  effect,  all 
curves  obtained  at  the  same  viscosity  ratio  should  reproduce 
themselves.  This  was  not  the  case.  Figure  11  shows  that 
higher  rates  favored  a  higher  water  fraction  flowing.  It 
was  noted  that  all  curves  extrapolated  to  fw  c  0  at  a  water 
saturation  between  11#  and  13#  which  is  a  reasonable  con¬ 
firmation  of  the  irreducible  water  saturation.  The  lowest 
rate  extrapolated  to  the  highest  water  saturation  in  line 
with  previous  comments  on  the  Jamin  effect  but  the  trend  was 
not  definitely  established  at  all  rates. 

The  pressure  profiles  provided  an  excellent  picture 
of  the  effect  of  fingering  on  the  saturation  distribution 
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within  the  core.  A  convenient  presentation  was  obtained  by 
plotting  the  pressure  drop  in  terms  of  the  inverse  of  the 
relative  injectivity.  The  influence  of  rate  and  viscosity 
on  the  pressure  distribution  at  a  stage  before  breakthrough 
is  shown  in  Figures  12  and  13  respectively.  Similar  curves 
are  shown  in  Figures  14  and  15  at  the  instant  of  breakthrough. 

It  may  be  noted  that  the  profiles  tended  to  flatten  out  at 
high  rates  and  viscosity  ratios,  indicating  viscous  fingering. 

The  effect  of  fingering  on  displacement  efficiency 
is  shown  in  another  way  by  plotting  the  injectivity  history 
with  rate  and  viscosity  ratio  as  the  parameters  in  Figures  1 6, 
III-5  and  17.  Referring  to  Figure  16  it  can  be  seen  that 
the  low  rate  injectivities  were  lower  than  the  high-rate  inject¬ 
ivities  prior  to  the  approach  of  breakthrough,  but  higher  than 
the  high  rate  injectivities  following  breakthrough.  The 
anomaly  of  the  lowest-rate  injectivity  being  greater  than  one, 
at  the  start  of  displacement,  required  an  explanation.  This 
was  believed  due  to  the  fact  that  as  soon  as  the  displacement 
commenced  the  measured  pressure  drop  A  P  differed  from  the 
initial  pressure  drop  A  by  the  capillary  pressure.  The 
value  used  for  the  initial  pressure  drop  was  always  that 
obtained  for  oil  flow  at  connate  water  saturation.  The  actual 
pressure  attributed  to  capillary  pressure  of  about  5.5  psi  was 
only  significant  at  the  very  lowest  rate  employed. 

In  an  effort  to  evaluate  a  possible  relation  between 
finger  growth  and  the  distance  travelled  by  the  front  along  the 
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core,  hypothetical  breakthrough  recoveries  were  calculated 
which  were  based  on  the  pressure  profiles.  Hypothetical 
breakthrough  recovery  is  that  which  would  be  obtained  if 
the  core  length  coincided  with  the  position  of  the  front  as 
revealed  by  the  pressure  profile.  Partial  pore  volumes  were 
calculated  assuming  a  homogeneous  pack.  The  result  of  this 
study  is  shown  in  Figure  18.  Although  there  is  considerable 
scatter  a  general  trend  can  be  observed  toward  higher  break¬ 
through  recovery  as  the  front  progressed.  It  appeared  that 
finger  growth  was  not  accelerating  but  rather  stabilizing. 
Because  maximum  breakthrough  recoveries  were  obtained  at 
the  full  length  of  the  model  it  was  concluded  that  still 
higher  recoveries  might  be  obtainable  in  a  core  longer  than 
eight  feet.  It  was  not  obvious  whether  rate  could  be  estab¬ 
lished  as  a  parameter  in  Figure  18  although  there  was  a 
general  trend  toward  stability  at  lower  rates. 

Subordinate  Production.  The  subordinate  production  data  was 
analyzed  to  determine  the  effect  of  fingering,  if  any,  on 
production  after  breakthrough. 

An  example  of  pressure  profiles  obtained  after 
breakthrough  is  shown  in  Figure  19.  It  was  noted  that  the 
first  two  profiles  obtained  following  breakthrough  still 
exhibited  a  frontal  appearance  normally  associated  with  a  high 
saturation  gradient. 

Following  a  procedure  similar  to  that  used  by  Loomis 
and  Crowell(36)  a  plot  was  made  of  f0  against  Qj_  on  semi  log 
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FIG.  19  EXAMPLE  OF  PRESSURE  PROFILES  AT  VARIOUS  CUMULATIVE  PV.'S  INJECTED  AFTER  BREAKTHROUGH. 
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paper  as  shown  in  Figure  20.  No  theoretical  basis  can  be 
given  to  justify  the  interpretation  of  a  linear  plot  indicated. 
This  is  a  constant  rate  of  decline  which  sometimes  is  actually 
observed  in  fields^2?) .  The  purpose  of  this  plot  was  to 
determine  a  cut  off  for  rejecting  data  which  might  be  influenced 
by  a  saturation  gradient.  The  existence  of  an  associated 
capillary  pressure  gradient  negates  the  application  of  the 
simplified  fractional  flow  formula  in  computing  relative 
permeability  ratios. 

Loomis  and  Crowell  were  concerned  with  the  existence 
of  the  "stabilized  zone"^^)^  This  is  the  region  in  which 
capillary  and  viscous  forces  are  postulated  to  be  in  balance. 

In  this  work  the  existence  of  a  stabilized  zone  would  not  be 
expected  except  possibly  at  the  very  lowest  rate.  The 
"frontal  zone"  of  Figure  20  is  interpreted  to  be  a  dispersed 
zone  of  high  water  saturation  somehow  related  to  fingering. 
Unfortunately  only  four  of  the  tests  were  carried  to  high 
enough  water  oil  ratios  to  permit  analysis  of  the  subordinate 
production  beyond  the  "frontal  zone". 

Theoretical  Recovery.  In  order  to  determine  whether  the 
recoveries  obtained  in  this  work  could  be  predicted  by 
Buckley -Leverett  theory  relative  permeability  ratios  were 
required.  The  use  of  the  unsteady  state  relative  permeability 
data  as  shown  in  Figure  5  was  not  entirely  justified  because 
those  data  were  derived  assuming  Buckley -Leverett  theory  did 
apply.  The  steady  state  relative  permeability  data  obtained 
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as  shown  in  Tables  II-l  to  II-4  of  Appendix  II  were  there¬ 
fore  obtained  . 

The  theoretical  recovery  as  a  function  of  viscosity 
ratio  is  shown  in  Figure  22.  The  recoveries  were  obtained 
from  the  fractional  flow  curves  constructed  in  Figure  21. 

Because  the  steady  state  relative  permeability  data 
did  not  extend  to  a  high  enough  saturation  range  it  was 
necessary  to  construct  a  composite  curve  based  on  all  data 
as  shown  in  Figure  23.  The  use  of  this  composite  relative 
permeability  ratio  curve  resulted  in  theoretical  recoveries 
less  than  actually  obtained  in  the  tests.  However  this  was 
not  considered  to  be  an  indication  of  the  inapplicability  of 
Buckley -Leverett  theory  to  this  case.  It  was  found  that 
theoretical  recoveries  were  extremely  sensitive  to  the 
position  and  shape  of  the  relative  permeability  ratio  curve 
in  this  saturation  range.  The  difference  between  actual  and 
theoretical  recoveries  can  therefore  be  explained  by  experi¬ 
mental  error  of  a  magnitude  which  is  suggested  by  the  spread 
of  data  points. 

To  further  evaluate  the  relative  permeability  data 
they  are  compared  with  those  of  Leverett  in  Figure  24.  The 

unsteady  state  relative  permeability  ratios  were  converted  to 
effective  permeabilities  by  the  method  of  Johnson,  Bossier 
and  Naumann^ 25) .  An  example  of  this  method  is  given  in 
Table  III-27.  Values  for  the  relative  permeability  to  water 
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for  this  work  are  slightly  lower  than  those  of  Leverett, 
particularly  the  unsteady  state  points.  The  values  of 
relative  permeability  to  oil  are  significantly  greater  than 
those  of  Leverett.  This  is  compatible  with  the  theory  of 
Odeh. 

Theoretical  recoveries  based  on  Leverett’ s 
relative  permeabilities  are  lower  than  those  obtained  from 
this  work.  This  is  shown  in  Table  3  which  compares  the 
recoveries  cited  in  the  literature  for  various  viscosity 
ratios  with  the  ones  derived  from  these  data.  Recoveries 
obtained  in  this  work  were  generally  higher  than  those  of 
other  workers.  Recoveries  obtained  by  Kyte  and  Rapoport 
were  significantly  higher  however  and  it  may  be  pertinent 
that  of  all  the  papers  referred  to  in  Table  3  theirs  was 
the  only  one  which  was  cognizant  of  the  possible  existence 
of  fingering.  On  the  other  hand  the  recoveries  shown  for 
Fatt  and  Klikoff  and  Leverett  were  derived  from  their  relative 
permeability  curves  and  therefore  should  not  be  influenced  by 
fingering. 
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CONCLUSIONS 

The  following  conclusions  were  made  on  the  basis 

of  the  test  results. 

1.  Breakthrough  recoveries  of  24.6$  and  29.8$  of  initial 
oil  in  place  were  obtained  at  viscosity  ratios  of  646 
and  292  respectively. 

2.  Actual  breakthrough  recoveries  obtained  were  higher 
than  those  predicted  by  Buckley -Leverett  theory.  The 
discrepancy  between  theoretical  and  observed  recoveries 
was  believed  to  be  due  to  the  extreme  sensitivity  of 
theoretical  recovery  to  small  changes  in  the  relative 
permeability  ratio  through  the  range  of  low  water 
saturation  involved. 

3.  There  was  some  evidence  that  higher  recoveries  would  be 
possible  with  a  longer  model. 

4.  Breakthrough  recovery  was  found  to  be  sensitive  to  the 
rate  of  displacement.  High  rates  had  a  detrimental  effect 
on  displacement  efficiency.  At  lower  rates  the  influence 
diminished.  This  rate  effect  was  attributed  to  viscous 
fingering  and  with  reference  to  the  theory  of  Chuoke  a 

C  value  of  25O  was  tentatively  assigned  to  the  model  to 
establish  the  point  of  onset  of  fingering. 

5.  Subordinate  production  was  divided  into  two  flow  regimes. 
The  first  regime  was  that  consisting  of  a  fairly  large  oil 
fraction  which  came  immediately  after  breakthrough.  The 
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second  was  that  of  a  smaller  oil  fraction  which  declined 
at  a  constant  rate.  Only  production  in  the  second  regime 
was  used  to  calculate  unsteady  state  relative  permeabilities 
because  of  the  high  saturation  gradients  associated 

with  the  first  regime. 

6.  The  sand  used  in  the  tests  was  found  to  be  preferentially 
water  wet  and  its  wettability  remained  unchanged  during 
the  tests . 

7.  The  high  viscosity  ratio  was  unfavorable  to  the  oil 
mobility.  This,  combined  with  the  saturation  change 
accompanying  the  passage  of  the  front  resulted  in  a 
sharp  loss  of  oil  mobility  behind  the  front. 

8.  The  'wet  sand  packing  procedure  resulted  in  a  model 
of  variable  permeability  throughout  its  length.  The 
permeability  increased  toward  the  base  of  the  pack. 

9.  The  variations  in  model  properties  between  the  two  packs 
did  not  appear  to  have  influenced  recoveries. 

10.  In  the  crude  oil  viscosity  range  of  170  cp.  to  1500  cp. 
original  brine  saturations  were  restored  at  the  conclusion 
of  a  displacement  by  injecting  as  little  as  one  pore 
volume  of  crude. 

11.  Using  this  resaturation  technique  as  opposed  to  cleaning, 
drying  and  resaturating,  in  conjunction  with  viscosity 
ratio  control  by  varying  temperature  reduced  the  time 
between  tests  from  about  one  week  to  one  day. 
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RECOMMENDATIONS 

Before  an  attempt  is  made  to  transfer  the  behavior 
of  this  model  to  other  systems,  some  verifying  tests  snould 
be  run.  An  attempt  should  be  made  to  further  examine  the 
effect  of  rate  on  the  oil  permeability  at  initial  connate 
water  saturation  in  the  absence  of  disturbing  factors  such 
as  temperature  change,  variable  Interfacial  tension,  and 
corrosion.  Due  to  the  scarcity  of  correlative  evidence  from 
other  unconsolidated  sand  models,  the  possible  influence  of 
those  variables  unique  to  this  model  cannot  be  discounted. 

These  factors  could  all  be  controlled  using  the 

same  model  in  the  following  way.  The  same  sand  could  be 

initially  treated  more  extensively  to  render  it  water  wet. 

The  use  of  an  additive  to  the  brine  such  as  sodium  tripoly- 
( 24) 

phosphate  would  help  reduce  the  corrosion  and  at  the 

(23) 

same  time  preserve  the  water  wetness  of  the  system  .  'y 
using  an  artificial  brine  and  a  blend  of  refined  mineral 
oils  the  interfacial  tension  could  be  controlled.  The  tests 
could  be  run  at  constant  temperature  at  a  viscosity  ratio 
determined  by  blending. 

The  selection  of  the  rate  of  injection  should  be 
based  on  the  rate  studies  of  this  thesis  with  reference  to 
the  particular  viscosity  ratio  employed. 

All  future  displacement  tests  should  be  carried  to 
a  water  to  oil  ratio  of  at  least  ten  to  one.  Subordinate 


■ 
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production  data  should  be  analyzed  in  an  attempt  to  determine 
if  the  point  of  separation  between  the  primary  and  subordinate 
phases  is  really  breakthrough,  or  some  time  later  in  the 
displacement  history. 

A  study  should  be  conducted  on  the  effect  of 
length  on  recovery  in  the  viscosity  range  under  study.  To 
increase  the  model  length  significantly,  a  length  of  thirty 
or  forty  feet  might  be  considered.  Pressure  taps  could  be 
placed  at  the  ends  only  but  provision  could  be  made  to  measure 
the  pressure  in  each  phase  for  low  rate  studies. 

End  faces  should  be  designed  to  permit  flow  both 
through  and  across  the  entire  end  of  the  core.  This  would 
appear  to  be  of  increasing  importance  as  the  model  diameter 
to  length  ratio  increases. 


V 
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PRACTICAL  APPLICABILITY 

The  results  of  this  work  indicate  that  conventional 
water  flooding  may  be  attractive  in  systems  similar  to  the 
model  in  the  viscosity  range  of  100  to  1000  centipoise  or 
even  higher.  However  the  unfavorable  mobility  ratio 
associated  with  the  subordinate  production  is  related  to  high 
injectivities  after  breakthrough.  Homogeneity  of  the 
reservoir  is  therefore  a  property  which  should  be  investigated 
thoroughly  prior  to  flooding.  If  stratification  exists,  any 
success  in  counteracting  bipassing  would  probably  yield  a 
more  dramatic  improvement  in  recovery  than  viscosity 


reduction  alone. 


NOMENCLATURE 


A  - 

cross  sectional  area  (cm2) 

C 

constant 

fo 

fraction  of  oil  flowing  at  any  point 

*w 

fraction  of  water  flowing  at  any  point 

I 

capillary  scaling  number 

X .  0 . 1 .  P  .  r 

■  initial  oil  in  place 

K 

absolute  permeability  ( darcies) (cm2) 

Kw,  Kq  - 

effective  permeability  to  water  and  oil  respectively 

(darcies) 

kV)J  ^0  ~ 

relative  permeabilities  to  water  and  oil  respectively 

L 

length  (cm) 

LvX,  _ 

scaling  coefficient  (£2£l£E) 

(  min  ) 

)un  _ 

average  distance  between  fingers  (cm) 

^cr  _ 

critical  wavelength  (cm) 

0 

porosity 

P.V.  _ 

pore  volume  (fraction) 

pc  - 

capillary  pressure  (atm) 

AP 

differential  pressure  (psig) 

(AP)l- 

differential  pressure  at  start  (psig) 

«1  - 

cumulative  injection  (P.V.) 

<*t  - 

total  flow  rate  (cc/sec) 

water  saturation  (fraction) 

swc  - 

connate  water  saturation  (fraction) 

:  -  ...  : 


■* 1 


o 

/U  w 
<T"  - 

V 

WOR  - 


oil  viscosity  (cp)(ps) 
water  viscosity  (cp)(ps)* 
surface  tension  (dynes/cm) 

total  flow  rate  per  unit  cross  sectional  area 
(cm/sec) 

water  oil  ratio 


* 


where  two  units  are  given  the  second  unit  applies  to 
finger  scaling  theory. 
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TABLE  1-1 


ANALYSIS  OF  RESERVOIR  WATER 


Appearance  of  Sample:  Clear  liquid 

Small  amount  of  rust-colored  sediment 


Percent  of 

Mg  per  liter  Calculated  Solids 


Cl  51, “67 

£8§3  118 

SO4  2 

OH  0 

Br  195 

I  12 

Ca  3,5^0 

Mg  1,681 

Na*(Calc)  26,2^3 

61.8 

0.0 

0.2 

0.0 

0.0 

0.2 

0.0 

4.3 

2.0 

31.5 

Total  Solids  83*258 

(Calc) 

100.0 

*  Alkali  metals  calculated 

as  Na . 

Density,  60°f 

1.0583 

pH 

6.2 

h2s 

None 

Resistivity,  ohm-meters, 
25°C 

0.101 

Ref.  Index,  25°C 

1.3^69 

Total  Solids 

evap.  at  110°C 
evap.  at  l80°C 

83,470 

83,090 

after  ignition 

81,830 

Remarks:  Fe-some 


1600 


7k  - 


S3SIOdllN3D  A1ISODSIA 


100  200  300  400  500  600  700  800  900  1000 

PRESSURE  PSIG 
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FIG.  1-2  VISCOSITY  OF  DEAD  LLOYDMINSTER 
CRUDE  AS  A  FUNCTION  OF  TEMPERATURE 


VISCOSITY  CENTIPOISES 
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STEADY  STATE  RELATIVE  PERMEABILITY  PRESSURE  PROFILES 
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TABLE  II -4 

STEADY  STATE  RELATIVE  PERMEABILITY  STUDY 
Sample  Calculation 

WOR  5:1;  25:5  cc/hr 

v  q  u^L 

K  =  Tap” 

8  x  30-5  x  14.7  /q  u\ 

11.43  x  3600  Iap/ 

=  0.087  (§^J 

where  q  is  in  cra^/hr 
u  is  in  cps 
A  P  is  in  psig 
K  is  in  darcies 


Bath 

Temperature  -  100 °F 

//*  -  0.810 

Mean 

Pressure  -  139.1 

yU  0  -  520  0PS 

A?  - 

108.2 

Kabs.  "  3.62  darcies 

Stabilized  Sw  -  0.288 

K  r  ,x.  .5,  ,x.  .5?0,  ^  2.09  darcies 

0  108.2 

k0  ~  -  0.578 

=  8,7  x  10 ~-2  x-  -25  x  Q>81'  -  0.0163  darcies 

^  108.2 

\r  =  °.<0^63  ^  0.0045 

3.62 

=  7.79  x  10 "3 
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